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Abstract 

We study stability of a circular orbit of a spinning test particle in a Kerr 

spacetime. We find that some of the circular orbits become unstable in the 

direction perpendicular to the equatorial plane, although the orbits are still 

stable in the radial direction. Then for the large spin case {S < 0(1)), the 

innermost stable circular orbit (ISCO) appears before the minimum of the 

effective potential in the equatorial plane disappears. This changes the radius 

of ISCO and then the frequency of the last circular orbit. 
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I. INTRODUCTION 



Coalescence of two compact objects, such as neutron star-neutron star, neutron star-black 
hole and/or black hole-black hole binaries, is one of the promising sources of gravitational 
waves which may be detected in the near future by a laser interferometric detector, such as 
US LIGO [|T|. If we detect a signal of gravitational waves emitted from those systems and 
compare it with theoretical templates, we may be able not only to determine a variety of 
astrophysical parameters of the sources, e.g., their orbital information, masses, and spins, 
but also to obtain more information about fundamental physics, e.g., the equation of state at 
high density In order to obtain such important information from the observed data, we 
have to make an exact template and then we need to know the exact motion of the binary. 
Hence the equations of motion in post-Newtonian framework have been studied by many 
authors We can see from the obtained equations of motion that spins of stars plays a very 
important role. For example, this effect induces a precession of the orbital plane, resulting in 
modulation of the gravitational wave forms 0|. Although such post-Newtonian approaches 
are definitely important, a fully general relativistic treatment must become necessary at 
some stage. Numerical relativity is one of most promising treatments to give a correct 
template, but it is still at a developing stage. Instead, a black hole perturbation technique 
has been well studied and provides us a good approximation. Then we expect that a general 
relativistic spin effect is also well understood by such a perturbation approach. 

Many studies about a relativistic spin effect using a spinning test particle around a black 
hole have been made since the basic equations were derived first by Papapetrou |^ and 
reformulated by Dixon 0. Bailey and Israel elaborated this system using the Lagrangian 
formalism 0. Corinaldesi and Papapetrou first discussed the motion of a spinning test 
particle in Schwarzschild spacetime §]. The Kerr or Kerr- Newman spacetime case was also 
analyzed by several authors 0. In |T^, the gravitational waves produced by a spinning 
particle falling into a Kerr black hole or moving circularly around it were discussed and 
the energy emission rate from those systems was calculated. We investigated more generic 
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motion of a spinning particle around a Schwarzschild black hole and pointed out that the 
spin effect can make some orbits chaotic |T^. 

When we discuss the gravitational waves from a binary system, one of the most important 
keys is the innermost stable circular orbit (ISCO). In the binary system, as gravitational 
waves are emitted, both energy and angular momentum of the system decrease, and so the 
orbital radius gradually decreases. When the angular momentum is reduced below some 
critical value, we cannot find any circular orbit. The binary system evolves into a more 
dynamical stage from the quasi-stationary stage. The last circular orbit is the ISCO. It 
may determine the observable frequency of the system. Using the effective potential of the 
spinning particle in the equatorial plane, the binding energy of the circular orbit is discussed 
in 0. They assumed that a stable circular orbit exists unless the minimum point of the 
effective potential in the equatorial plane disappears, finding the largest binding energy. 

However, before the system reaches such a last orbit, if there exists any unstable mode, 
it will change the ISCO, and therefore the last circular radius, binding energy and frequency, 
which may be important in gravitational wave astronomy. In fact, we found some circular 
orbits of a spinning particle in Schwarzschild spacetime become unstable in the direction 
perpendicular to the equatorial plane []T2[ . 



In this paper, extending our previous analysis ||T2|, we study an instability of the spinning 



test particle in a Kerr spacetime. The reason for a transition from the quasi-periodic stage 
to dynamical stage, we have so far three effects. Most important one is general relativistic 
strong gravity, i.e. when the orbital separation of the binary becomes small as r ~ 6M , the 
centrifugal force cannot be balanced with the gravitational force, so that the star cannot stop 
from changing its circular orbit to a plunging one. The second effect is the hydrodynamical 
interaction between two stars, i.e. at the moment when the surfaces of the stars come into 
contact, we expect that such a transition will occur due to the direct interaction between 
stars. The third one is an effect due to a property of the star such as a spin or a multipole 



moment. For example, in |]T3[, they show that the instability can occur due to a tidal 



force. When the orbital separation of the stars becomes small, each star of the binary is 



significantly deformed by a tidal force. As a result, a circular orbit of the stars becomes 
unstable and the ISCO appears earlier. In this analysis, however, the motion of the stars 
is restricted into the equatorial plane and the stars plunge to each other when the motion 
becomes unstable in the radial direction. Here, we look for a new type of instability due to 
an effect belongs to the third category, i.e. the instability in the direction perpendicular to 
the equatorial plane. We will show that such an instability really occurs when the spin of 
a test particle is high. A spinning particle leaves the equatorial plane and falls into black 
hole. Since it is important to determine the ISCO for the gravitational wave astronomy, we 
believe that it is also worthwhile to analyze such a new instability of a spinning particle in 
Kerr spacetime. 

This paper is organized as follows. After a brief review of the basic equations in section 
2, we analyze stability of a circular orbit of a spinning test particle in section 3, using a 
linear perturbation method. We show changes of the radius of the ISCO, of the binding 
energy and of the frequency due to such an instability. Summary and some remarks follow 
in section 4. Throughout this paper we use units c = G = 1. Our notation including the 



signature of the metric follows MTW [P 



II. BASIC EQUATIONS FOR A SPINNING TEST PARTICLE 

A. Pole-dipole Approximation 

The equations of motion of a spinning test particle in a relativistic spacetime were first 
derived by Papapetrou and then reformulated by Dixon . These are a set of equations: 

^ = ^^ (2.1) 

as 

Dn^ 1 

= -l^R'u,yS^% (2.2) 

^^=pt^v'^ -p'^yt^^ (2.3) 
ds 

where s,v^,p^ and S^^ are an affine parameter of the orbit = x^(s), the 4-velocity of 

the particle, the momentum, and the spin tensor, respectively. This is called the pole-dipole 



approximation, where the multipole moments of the particle higher than mass monopole and 
spin dipole are ignored. We need a supplementary condition which gives a relation between 
v'^ and p^, because is no longer parallel to v'^. The consistent choice of the center of mass 
provides such a condition 0, 

Pf.S^" = 0. (2.4) 



Using ( |2.4| ) we find the relation between v'^ and as 

1 



= N 
where 



(2.5) 



5 = 1 + -l^R^f^.^sS'^f'S^', (2.6) 

and is a normalization constant, which is fixed by a choice of the affine parameter s. 
= p'^/ /X is a normalized momentum, where the mass of the particle /x is defined by 

/i' = -Pupr (2.7) 

It may be sometimes more convenient or more intuitive to describe the basic equations 
by use of a spin vector S^, which is defined by 

S, = ~\t,,p.u''SP\ (2.8) 

where e^'^''" is the Levi-Civita tensor. The basic equations are now 

Dv^ 1 

^ -R^'^^^XS^P", (2.10) 



where 



ds jj 
n Qf^ 1 

— = -p'^R^.^^S-'v'S^p'^, (2.11) 



Hi/pa — 2"^^"^ ^^aPpa- (2.12) 
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Eq. (pl^) with the definition ( pl8|) reads 

PpS'' = 0, (2.13) 
which gives the relation between and p^, 

= u'^ + (2.14) 

where 

by fixing the affine parameter s using the condition v^u^ = —6. This choice makes the 
perturbation equations simpler as we will show later. 

For a particle motion in a Kerr spacetime, we find several conserved quantities. Re- 
gardless of the symmetry of the background spacetime, it is easy to show that /i and the 
magnitude of spin S, defined by 

S' = S^S>', (2.16) 

are constants of motion. When the spacetime possesses some symmetry described by a 
Killing vector C,^, 

C ^ - l^.-^uS^" (2.17) 

is conserved . For the spacetime which has both axial and timelike Killing vectors such as 
Kerr spacetime, we have two conserved quantities, i.e. the energy E and the z component 
of the total angular momentum of a spinning particle J^. For Schwarzschild spacetime, the 
"x" and "y" components of the total angular momentum of a particle are also conserved 
because of spherical symmetry. 

B. Effective potential of a spinning particle on the equatorial plane 



In order to discuss the motion of a test particle in a black hole spacetime, we usually 
introduce an effective potential. However, if the test particle has a spin, it is not so easy to 
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find an effective potential because of an additional dynamical freedom of the spin direction. 
In our previous paper, we defined an "effective potential" for a spinning test particle in 
Schwarzschild spacetime, i.e. V(±){r,9; J,S) given by Eqs.(2.31) and (2.32) in [Q. Setting 
the total angular momentum as to point in the z direction, i.e. {Jx,Jy,Jz) = (0,0, J), 
those equations are obtained from the condition oi = = in ( |2.4| ) and ([^.T]). We use 
quotation marks because it is not a real effective potential but plays a similar role. 
The particle with energy E can move within the contour curve defined by 

E = V^^){r,9-J,S). (2.18) 

This "effective potential" has two parameters, J and S, whose values determine the topology 
of its contours. From the contours of the "effective potential" (Fig.2 in |jl2[); we can easily 
see the stability of a bound orbit. 

However, since we are interest in a Kerr background spacetime, whose metric is given as 

ds^ = - f 1 - ^) dt' + ^dr^ + ^dO^ + ^ sin^ ed<p' - sin^ edtdct>, (2.19) 

where, 

E = + cos^ ^, A = - 2Mr + ^ = (r^ + a^)^ - a^A sin^ ^, (2.20) 

it is hard to define an effective potential even in the above naive sense. This is because we 
cannot find new conserved quantities in a Kerr background spacetime, which correspond to 
Jxi Jy for Schwarzschild case or Carter's constant for a spinless particle |]14|, in addition to 



E and J^. However, if we restrict the particle motion in the equatorial plane, an effective 
potential for the radial motion is found |0. The direction of spin must be parallel to the 
rotational axis. In fact, setting 6 = 7r/2, p"^ = 0, S'^ = = = and = —5, we find 
that Eq. ( p.7[ ) is reduced to, 

(pi)2 = A(E-[/(+))(i5; -[/(_)), (2.21) 

and 



[/(±)(r; J„ 5, a) = XJ, ±JiX^- Y).P, - Z, 



(2.22) 
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Note that the suffices of the momentum and the spin vector denote those tetrad components. 
The tetrad frame has been defined as 

—,0,0, -asm e\l — 

(2.23) 





ej= (0,0,v^,0 



sin ^,0,0, 



sin 9 



where = (e^, e°, e^, e^) for a = ~ 3. We use of S including the direction of spin as well 
as the magnitude. Then aS > means the spin of the particle is parallel to that of the 
black hole and aS < means anti-parallel. 

We regard f/(±)(r; J^, 5, a) as an effective potential of the particle on the equatorial 
plane. When a vanishes, [/(+) is reduced to V{r,7i/2; J, S) in the Schwarzschild case [p!2[ . 
The particle with energy E can move only in the region of -E > (or E < on 

the equatorial plane. The typical shape of is shown in Fig.l. As for the orbit in 
the equatorial plane, we usually discuss it by use of this effective potential. In the case of 
a spinless particle, it is enough because the angular momentum conservation prevent the 
particle from leaving off the equatorial plane. If the angular momentum is larger than some 
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critical value, we always find a stable circular orbit. Since a circular orbit does not exist 
below the critical angular momentum, the particle eventually enters a dynamical stage from 
a quasi-periodic stage because the emitted gravitational waves bring the angular momentum 
of the system away. This critical value gives the ISCO. If a test particle has a spin, however, 
the orbital angular momentum needs not be conserved. The particle may leave off the 
equatorial plane. Hence the analysis by use of the above effective potential may not give the 
ISCO. We need a more detailed analysis to see whether the circular orbit on the equatorial 
plane is really stable or not. 
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Fig. 1: The effective potential [/(j.) on the equatorial plane for = 4/xM, S = 1/xM 
and a = 0.5M. (a) The potential C/(+) has two extremal points, (b) The potential 
has no extremum. A particle with positive energy is not bounded. 

III. STABILITY OF A CIRCULAR ORBIT OF A SPINNING PARTICLE 



The linear perturbation method is used to analyze the stability of a circular orbit of a 
spinning particle. In the tetrad frame, the equations of the motion of a spinning particle are 
rewritten as 



ds 
dp"" 

ds 

dS^ 

ds 



a„,c nb 



S' + -p-S'R*,^^^v^S'p^ 



(3.1) 
(3.2) 
(3.3) 



where Uijk are the Ricci rotation coefficients defined as 

^^ijk = eic'^ek^.^. (3.4) 

We assume a circular orbit in the equatorial plane as an unperturbed orbit, i.e. r(o) = 
ro = const., 6'(o) = it/2, pj^-^ = pf^^ = 0, S'J,) = Sfg) = S'gj) = 0, S'g^^ = -S and the 
time derivatives of these variables vanish. Here, we use the suffix (0) for the unperturbed 
variables. From Eg. ( p.l4| ), we find 

and f(Q) = f^Q) = 0. The conserved quantities in our system are the energy and the z 
component of total angular momentum, which are now given as 



The way to determine the unperturbed variables is as follows. First, the Kerr parameter a 
and the magnitude of the spin S are fixed. Next, we give the radius of a circular orbit tq. 
Then, from the condition of potential extremum at tq, i.e. 

dVi+){r) 
dr 

r=ro 

Jz is determined. The energy of the particle is obtained hj E = V(+)(ro). Finally, from 



0, (3.9) 



Eqs.( |3.7| ) and (p^), p^Q^ and p^Q^ are determined. We denote these unperturbed variables as 
^(0)' "^(0)' ^'(0) ^{0)- Of course, if no real solution is found for given parameters a and 
iS, a circular orbit does not exist at tq. 

Let 6x^, 5v\ and 5S'^ be the perturbations around the unperturbed variables. We 
shall introduce a tetrad description for as = e^^5x^. From Eq. (|2.14|) , 

5v' = Xpx^ + Yj6p^ + Zi5S^ , (3.10) 
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where, 



Inserting 



1 (9 

u__*p* r ql cm n 
— -"-(0) «mn'-'(0)^(0)^'(0)> 



1 1 

""i + -"(0) imn'^(0)'^(0)fi ) 



3*^^0) ^mnl"i^(0) 



(^.'^(^)+^(o)<^r)P(o)- 



(3.11) 
(3.12) 
(3.13) 



•^(0) + i^^; , 
^(0) + -^^'^ 



P(o) + ¥ 



(3.14) 



into Eqs. (|3.1|) - (|3.3|) and neglecting higher terms than the first order of perturbed variables, 
we find the hnear perturbation equations as 



where 





5x* 




d 

ds 


dp" 











A^, B^, ^ ^ 
D^^ E^ F) 



r^i ZJi T'i 



(3.15) 



5* 



J' 



d 



* ^k 



'^mk%)Sj + ^mk'y'jPio) + n 



(O)Zfc 



J' 



^ ilk i I k i I 



(3.16) 



7§"-^{0)a6cd'S'('o)'i^('o)'S'(o)P^O)'^i + ^(0)«fe*^j'S'('Q) + — p*S'(o)in' j-, 

■ I ^ ■ b : I k ^ ■ I 

^(0)lj *'^(0) + J^P\o)R*0)jabc'^{0)S(o)P{0) + ^{0)lk''^ jS(o) + ^P(0)'S'(0)/^ j, 
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with 



ff^i ^ I ^/J ^ p* i „,l I p* i vl \ om n 

^ i ~ ~ 1 ^i^^(O) «mn^(0) + ^(0) Imn^ j '^(0)^(0)' 

= —Rlo)\mn(X^jP"o) + ^{0)^7)'^{0)' (3-17) 

A'' 



Eq.( |3.15|) is nothing but the equation of geodesic deviation for the case of a spinless particle. 



The eigenvalues of the matrix in Eq. ( ^.151) determine the stability of a test particle. Inserting 



the unperturbed variables into the matrix in Eq. (|3.15|) , the matrix becomes rather simple. 
In Appendix, we show the explicit form the matrix and its components and discuss the 
structure of the matrix. The eigenequation forms 

X%X^ - Ar){X' - AeiX^ + Ae2) = 0, (3.18) 

where A^, Aei and Ag2 are defined in Appendix. First, to see the meaning of this equation, 
we consider the simple case . When S = and a = 0, this equation is reduced to 



I.e. 



^ Mfx\ro-6M) 
rl{ro-3M) ' 

Aei = (3.20) 





r, 



Ae2 = 0, 



where Lz is the z component of orbital angular momentum, and it is given for an unperturbed 
circular orbit as 

^^-ro-3M- ^^-^^^ 

This equation shows that any circular orbit does not exist for tq < 3M. The positivity of the 
second parenthesis of (|3.19|) guarantees the stability in the 6 direction, i.e. in the direction 
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perpendicular to the equatorial plane. Then, the orbit of a spinless particle in Schwarzschild 
spacetime is stable against perturbations in the 9 direction. The first parenthesis in ( |3.19| ) 
shows the stability in the r direction. For tq < 6M, the eigenvalues are real, while for 
r > 6M, those are pure imaginary. This means that the equilibrium point (circular orbit) 
with ro < 6M is unstable, while it is stable for the case of rg > 6M. We find the critical 
radius r = 6M when Lz = \/l2fiM. This corresponds to the ISCO of a spinless particle in 
Schwarzschild spacetime. 

In the general case, we can analyze stability of the circular orbit in the r direction by A,, 
and in the 6 direction by Agi and A02 (see also Appendix). If A^. > 0, the circular orbit is 
unstable against a radial perturbation, while if a real part of any solutions of the equation 

- AeiX^ + Ae2 = 0, (3.22) 

is positive, it is unstable against perturbations perpendicular to the equatorial plane. 

In Fig. 2, we show the unstable parameter ranges by shaded regions for the case of a = 
0, 0AM and 0.8M. We find the unstable parameter range is usually narrow, because we 
have a constraint for a particle spin as S/fiM < 0(1) (see discussion in [0). Although the 



generic features of the parameter range is similar, there is some small difference depending 
on a. As we can see from this figure, the unstable region in the a = 0AM case is larger 
than that in the a = case. The value of S at the bottom edge of this region is ~ 0.9/iM. 
However, in the a = 0.8M case, the unstable region becomes narrow again. The lowest value 
of S for this region is about 0.94/iM. 

Fig. 3 shows the radius of the ISCO for several values of a. The radius of the ISCO 
decreases as S increases or as a increases, which means that the spin effect plays the role 
of a repulsive force for aS > basically. Therefore, the particle can approach the horizon 
of the black hole enduring the gravity. But, if S is larger than the critical value, the radius 
of the ISCO increases. This shows the occurence of the instability of the motion in the 9 
direction. In Fig.4, we show the energy of the particle and the orbital frequency of the ISCO. 
In 1^, the energy the ISCO for same system was analyzed, but only the instability in the 
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radial direction was taken into account. When the instabihty in the 9 direction is included, 
the energy and frequency of the ISCO behaves in different manner beyond the critical value. 



a = a = 0.4 [M] a = 0.8 [M] 




3.6 3.7 3 3.1 2.05 2.15 

J, [nM] [nM] J, [nM] 

(a) (b) (c) 

Fig. 2: The parameter regions for which the circular orbit becomes unstable in the 
direction perpendicular to the equatorial plane. The minimal values of the spin for 
each case are (a) S = 0.98fiM (b) S = 0.9^M and (c) S = OM^iM 




-1.5 -1 -0.5 0.5 1 1.5 
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Fig. 3: The radius of the innermost stable circular orbit (ISCO). Above the critical 
value of S, the radius of the ISCO increases. 
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Fig. 4: (a) The particle energy E of the ISCO. (b) The orbital frequency of the ISCO 
for which the mass of the black hole is assumed to be IOMq 

IV. DISCUSSIONS AND SUMMARY 



In our previous paper [T^, we discussed the motion of a spinning test particle in a 



Schwarzschild spacetime by use of an "effective potential". The "effective potential" is 
classified into four types depending on its topology of contour (Figs. 2 and 3 in [|12[)- From 



Fig. 2 in |T2|, we can see easily the stability of the circular orbit of the spinning particle in 
Schwarzschild background spacetime. The type(Bl) has one saddle point and one minimal 
point on the equatorial plane. The saddle point is maximal in the r direction and minimal 
in the 6 direction. At the minimal point, it is minimal in both directions. For the type(B2) 
potential, there are two saddle points apart from the equatorial plane and a maximal and 
minimal points also exist on the equatorial plane. In the type(U2) potential, we find one 
unstable point in both directions and one saddle point. At this saddle point, we find that a 
particle is stable in the r direction but unstable in the 6 direction. This is exactly the same 
as what we found in this paper as new type of instability. The stability of the circular orbit 
on the equatorial plane is summarized in Table 1. 

For a Kerr background spacetime, however, we cannot define the effective potential for a 
spinning particle like Schwarzschild case. Thereupon, instead of the topology of the contour, 
we divided the parameter space by the stability of the circular orbit. The result for a = 
case is shown in Fig. 6. The name of classified types in this figure are same as those of the 
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"effective potential" in in which the circular orbit shows the same stability type as the 
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present analysis. Since the region of > 3.5/iM and S > coincides with Fig. 3 in |]T2| 
believe that the present classification is the same as the previous one by use of the "effective 
potential". Fig. 7 is the classification for the Kerr background spacetime. We expected that 
the parameter region for type(U2) potential for Kerr background spacetime would be larger 
than that for Schwarzschild spacetime because of the lower symmetry and the existence of 
the spin-spin interaction. However, as we have seen in Fig. 2, the parameter region is not 
enlarged abruptly as a increases. Rather, for large a, this region gets narrow. 



Since our system is not integrable as shown in |T2|, it may be interesting to mention 
chaos of a spinning test particle in Kerr background spacetime. In our previous paper, we 
conclude that the motion of the particle can be chaotic in the type(B2) potential. As we see 
the type(B2) region becomes larger as a increases. Therefore, we expect that the motion 
of the particle in Kerr spacetime may become chaotic more easily than in Schwarzschild 
spacetime. This is consistent with the fact that Kerr spacetime possesses lower symmetry 
than Schwarzschild spacetime and that number of the constants of motion is smaller than 
that in Schwarzschild case. As for a example, we show in Fig. 8 the Poincare map of the 
chaotic orbit for the a = 0.8M case constructed with the same rule in |T^. The parameters 



of the orbit are chosen to Jz = 4yuM and S = 0.3/iM, which are in the type(B2) region. For 
Schwarzschild case, the magnitude of the spin at the bottom edge of the type(B2) region 
with Jz = 4:^M, i.e. the lowest value of the magnitude spin for which the chaos occurs, 
is about 0.635/iM (Fig. 7 in |jl2[|), which is mush higher than that of the above example 



{S = 0.3/im). This result supports our expectation [jT5 . 

In this paper, we have analyzed the stability of a circular orbit of a spinning test particle 
in a Kerr background spacetime using linear perturbation method. The circular orbit can 
be unstable not only in the r direction but also in the 6 direction due to the spin. With 
the parameters within the shaded region in Fig. 2, this instability occurs. In Fig.5, a typical 
unstable orbit with these parameters and its fate are shown. We have chosen the parameters 
for the initially circular orbit of the particle as 5 = l^M, E = 0.86199950863/i, Jz = 
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3.06620962243//M, vq = 3.05M and 9 = n/2. The Kerr parameter is a = 0AM. As for a 
perturbation, the direction of the spin is set to be shghtly different from the perpendicular to 
the equatorial plane, i.e. SS^ = 3.672835 x 10~^/xM. Because this circular orbit is unstable 
in the 9 direction, the particle leaves the equatorial plane and falls into the black hole 
eventually. Then this type of instability may change the radius of the ISCO. In Fig.3, the 
dependence the radius of the ISCO upon the magnitude of the spin is shown. In Fig. 4, the 
energy of the particle on the ISCO and the orbital frequency of the ISCO are also depicted. 

If this new type of instability occurs in a real astrophysical system such as a binary 
system, it changes the radius of the ISCO. This is very important for observation, especially, 
gravitational astronomy as mentioned in Introduction. So, we may restrict the parameters 
of the system, for example the lower limit of the spin, from the observation of gravitational 
waves. And the chaos is also important factor for the gravitational wave astronomy, because 
if the chaos occurs in the system, it must affect the gravitational wave emitted from the 
system. However, the test particle model should not be adopted as a model of the last 
stage of a binary system. In fact, in our analysis, the radius of the ISCO gets very small. 
For such small distances to the horizon, instabilities due to another effects mentioned in 
Introduction must be taken into account. Therefore, we need further investigation to check 
the importance of this type of instability in a more realistic binary system. 
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Fig. 5: The orbit of a spinning particle in type(U2) potential in Kerr background 
spacetime with a = 0.4M. The particle is at the minimum point of the "effective 
potential" in the r direction on the equatorial plane and the direction of the spin vector 
is slightly deviated from the direction perpendicular to the equatorial plane initially. 
The parameters of the orbit are = 3.06620962243;uM, E = 0.86199950863;U, S = 
l.OfiM and tq = 3.05M. The initial perturbation is 6S^ = 3.672835 x 10"^^- 
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Type(B2) 
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Fig. 6: The classification of the "efi'ective potential" of the particle for the a = case 
in the Jz-S plane. The region > 3.5fiM and 5 > is the same as Fig. 3 in 




Fig. 7: The classification of the effective potential of the particle on the equatorial 
plane for a Kerr black hole. 
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0.015 




6.3 6.4 6.5 6.6 6.7 

r [M] 

Fig. 8: The Poincare map for the orbit in a Kerr spacetime, a = 0.8M, with parameters 
which belong to the type(B2). All orbits has = 4.0 fxM, S = 0.3/iM and E = 
0.958155/x. We set p'' = and S'^'^/S^'^ = tan(85°) = 11.4301 initially. And the initial 
position is ro = 6.5, 6.592, 6.628 and 6.68M. We find the motion is chaotic. 





inner extremal point 


outer extremal point 




r direction 


6 direction 


r direction 


6 direction 


Type(Bl) 


Unstable 


Stable 


Stable 


Stable 


Type(B2) 


Unstable 


Unstable 


Stable 


Stable 


Type(Ul) 


no extremal point 


Type(U2) 


Unstable 


Unstable 


Stable 


Unstable 



TABLE. 1: The stability of a spinning test particle on the equatorial plane. 
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APPENDIX A: THE PERTURBATION MATRIX AND ITS EIGENEQUATION 



Here, we present the exphcit form of the perturbation matrix in Eg . (13 .151) and its non- 
vanishing components. 
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IZE Pio) + Z2^(0) + -:^^(o)> 



3MS 



(^fo)P?0) + ^{0)Ph) Pfo) - -2^?0) 



A 3 



A prime denotes d/dr. It is easily shown that the eigenequation of the matrix 
decomposed as follows. 



1|) can 



-A B\ 

-A 

-A E\ 
D\ E\ E\ -A 



-A D\ F\ 
-A G\ I\ 
-A G^^ I\ 

B\ C\ C\ -A 



H\ I\ I\ 



-A 



(A2) 



Three zero eigenvalues comes from two zero columns and a zero line in the matrix, which 
are related to the facts that t and are the cyclic coordinates of the system and that the 
magnitude of the spin is conserved, respectively. The second part consists of the compo- 
nents related to 5x^ , 6p^, 6p^ and 6p^. Because of the conservation of the energy and the 
angular momentum of the particle, two eigenvalues of this part are zero. The rest two eigen- 
values describe the stability in r direction. The components of the last part are related to 
Sp"^, 55*°, 6S^, 6x^ and 6S^. The eigenequation of this part, except for one zero eigenvalue 
which comes from the condition (|2.13|) , turns out to be a quadratic equation of A^ and 
determines the stability in the 6 direction. Finally, we find the eigenvalue equation as. 



A^(A2-A,,)(A^-AeiA' + A 



92 



0, 



(A3) 



where 



A, = {B\D\ + E\E\ + E\E\), 
Aei = {B\D\ + C\G\ + G\G\ + F\H\ + + 

A _ (r"^ p2 ^0 Tjl I ^2 172 ^3 ttI ,^2 7-)2 ttI tO 
A92 - gi* 2-n 2 + ^ 3-^ 1^ 2-n 2 ~ ^ 0^ 2-" 2-' 1 

r2 p2 ^0 t1 I r2 n2 rO rl , ^^2 ^3 jO jl 



(A4) 
(A5) 
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— £J 1*-^ 2-' 3 ~ '-^ 0*-^ 2-' l-* 3 ~ ^ 3-^ 2-" 2-' 1 

-C^3G°2/V\ + B\D\I\I\ + C\G\I\I\). (A6) 

To analyze a stability of a circular orbit, we have to know for the stability in the r 
direction, while A^i and A.02 for that in the direction. 
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